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Abstract 


Application of superconductines elect ricalemach merry ane 
ship propulsion yields light weight and small volume propul- 
sion plants with increased system efficiency, reliability 
and lower cost. Such propulsion systems might be the only 
feasible answer for high performance ship types. 


A desaen method for two=polée synchronous machine wig 
Superconducting field Winding is developed which vill Survive 
even under the worse fault conditions. Iron shielded machines 
Zach the effects of higher harmonics and partial winding 
angles are considered. The method is then adapted to a 
Computer program for ease of obtaining minimum volume апа 
ПІРІ machines. The results are compared with 20,000 HP 
conventional propulsion generator. it is shown that great 
weight and volume Savings are possible. 
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NOMENCLATURE 


Subseripts 

a armature winding 

f field winding 

р radial 

Ө azimuthal 

1 inside (for field quantities--at a radius less 
than inner radius) 

1 inner (when used with a dimension) 

O outside (when used with a field quantity) 

O outer (when used with a dimension) 

5 field quantities, generated by a sheet current 

a phase a 

b phase b 

E phase c 

d direct axis-armature 

q quadrature axis-armature 

sd direct axis-damper shield 

sq quadrature axis-damper shield 

Symbols 

B ПЕРИ ер ге ах density 

ЕВА saturation flux density of shield material 

C mean-squared field experienced by armature 
сошаиссог 5 

а, diameter of primary conductor 

den dissipation density of shield material (watts/1b) 


i 





open circuit voltage (RMS) 

voltage behind subtransient reactance 
voltage behind subtransient reactance (per unit) 
magnetic field density 

current 

rated armature current (RMS) 

current density 

rated armature-current density 

rated field-current density 

shield structure thickness 

primary shield thickness 
straight-section length 

Local machine "ση 


елисе Гепреиштокг calculating tar mature 
self inductance 


effective length for calculating field-winding 
self-inductance 


effective length for calculating field-to- 
armature mutual inductance 


effective length for calculating damper self- 
inductance 


effective length for calculating damper-to- 
armature mutual inductance 


Ν.Ε lengch tor calculating tield=to= 
damper mutual inductance 


ПО of armature conductors 


Su раје length тор calculating eddy-current 
losses 


effective length for calculating shield losses 


self inductance, armature phase 
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self inductance, iela 

mutual inductanee, phase-to-phase 

mutual inductance, armature-to-damper 

mutual inductance, field-to-damper 

magnitude of mutual inductance, field-to-armature 
mutual inductance, field-to-phase-a 

magnitude of mutual inductance, damper-to-phase-a 
number of turns ir armature winding 

number of turns іп field winding 

number of pole pairs 

rated volt-amperes 

conduction losses 

eddy-current losses 

shield dissipation 

inner radius, shield 

outer radius, M 

outer radius, damper shield 

Paced Cermna l Vor are 

synchronous impedance (in ohms) 


per-unit synchronous impedance, normalized to 
internal voltage, Er 


per-unit synchronous impedance, with V 
base voltage 


+ 25 


per-unit transient reactance (direct axis) 
per-unit subrtansient reactance 
БО О 


Βρι/ Βρο 








SKIN ας σα 

angular displacement 
fleld-winding included angle 
Pet 

armature winding angle 

Ре на 


magnitude of maximum radial stress under three 
phase fault from load 


radial stresses 

Sheer stresses 

density 

conductivity ОЕ НЕ ПСО 
conductivity of damper-shield material 
rotor-axis displacement from axis of phase-a 
angular frequency, also used for derlecittons 


electrical angular frequency = pw 
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INTRODUCTION 


Advantages of electrical ship propulsion systems Were 
realized years ago. Some of them are: arrangement flexib- 
ility, ease of automation and possible integration to ship 
service power. Major drawbacks were the relatively high 
eee and weight of system components. 

Over the past ten years or so many studies have been 
done on the possible applications of superconductors to 
electrical power equipment, with most effort devoted to 
nuansformers with deereasing attention given to transmis Tons 
ae machines, and de machines respectively. The development 
ee nigh-field superconducting wire, suitable forf bricat ion 
асо windings for electric machines, high vacuum seals and 
rotating dewar technology, made applications of supercon- 
итоге possible [ог кога есер? ы шав о 

Among the many different types of superconductors., ще 
Best useful One is type II, beeause of 1:65. ари πο TOR ке 
en current density in a high magnetic ela 
first major consequence of this upon machine design is that 
we may break free from the magnetic limitations of iron. 
INE mination of iron will reduce the machine weight and 
provide more armature conductors for a given space, thus 
mMer easing the power rating per unit volume. r sr ΠΠ 
density at the armature conductors is increased by a factor 
of 3 over that in a conventional machine, the generated emf 
Mer unit length of conductor is thus increased 3 Times, if 
шие armature diameter and the rotor speed are left unchanged. 
шизтегоге, Гот the same output power and current, the 
machine length can be reduced by the factor 3, with a 
corresponding reduction in machine volume and weight. 

Application of superconductors to rotating machines is 
not without problems. First of all, they will operate only 
at extremely low temperatures. Secondly, perhaps the most 


important one is that time varying magnetic fields produce 
iat 





losses which are prohibitive at power frequencies. These 
losses are the limiting factor for the application of super- 
conductors to field windings: Therefore те сотто Ее 
on a superconducting machine must be a superconducting 
stationary or rotating field winding and a normally conduct— 
ing ac armature winding, which may also be rotating or 
stationary. Configuration for rotating armature has inherent 
problem of transmitting power across sliding contacts: 

With the advent of rotatable liquid helium transfer couplingl, 
a rotating superconducting field winding became practical, 
thus eliminating the necesity for rotating the armature. 
Different designs and configurations have been studied in 


pecu 


aene rent countries, among them synchronous machines 
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de machines”, homopolar а m double armature ma- 
chines”. and the toroidal reciprocating alternator“. Some 
of them are already built or under development. 

Application of superconducting electric machinery to 
ship propulsion systems has been studied over approximately 


the past seven years /»9,9,10,11 


A superconducting propul- 
Био ο shown 1n Fisure 1. The reduir- di Sertion 
pertormanee characteristics from a propulsion plant wirl 
depend mainly on the ship type being considered. Such 
Aharacteristies are tied with the options, control and type 
of electric machinery. For example, bulk carriers which 
run at their full speed 85% of their operating time will 
not require speed control as much as a naval combatant 
vessel, and the option box may be omitted and a simple 
synchronous-synchronous drive in which the propeller speed 
can be controlled by varying prime mover speed, may be 
used. Destroyers, which operate over theirfull speed range, 
will require more sophisticated controls and the option box 
may be filled with a cyloconverter if a synehronous system 
is used or de machinery might be more appropriate. It 
should be pointed out that the choice of thruster and prime 


mover will aiso affect the electric machinery. 


l2 | 
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Recently some studies have been done on Conventions 


12,13,14. Sixteen 


electric propulsion of bulk carriers 
tankers, container ships and LNG carriers are being built 

ог аге under design. Conventional electric propulsion 
systems have been selected due to the fact that flexibility 
i Oropulsion systems may Lead to possible Savings in τας 
volume of the plant. Gas turbine generator sets can be 
placed on the main deck, reducing the inlet and exhaust 
mecing volume and flow pressure drops, thus improving gas 
turbine performance. Added to Chese advantages, ат а super- 
conducting drive is used even more savings are possible. 

It can be shown that? the reduction of рест л септе 
weight from 30 to 10 1bs/SHP will permit an increase in 
Maximum speed from 35 to 50 knots with no sacrifice in 
Payload or endurance for a destroyer type ship. Introduction 
of gas turbines has made specific weights somewhat lower, 
but introduced new probiem areas. Newly developed ship 
types such as hydrofoils, SES and catamarans which require 
unusual propulsion plant arrangements, suffer greatly with 
conventional propulsion plants. With increased power levels 
required by high performance ship types, superconducting 
electric propulsion plants may be the only feasible answer. 
The purpose of this thesis is to obtain a design method for 
Synchronous machines which can be used to realize possible 


advantages. 
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Chapter 1 
DEVELOPMENT OF A MODEL 


Model Configuration 


The geometrical model used for this analysis is 
melustrated im Preure 2 and Figure 3, with approximate 
account taken of the effects of end turns as indicated by 
measurements made on the first MIT experimental све i 
We akial length of one Set Of armature eng turns 15S ΕΤ 
to be equal to the average armature winding diameter in the 
active section divided by the number of pole pairs, and the 
radial chickness of the ena turns is assumed equal to twice 
τας radial thickness of the armature winding in the active 
шото nese relatively large end turnare assumed 
because there is more copper in the active region and more 
volume is needed in which to cross the phase conductors 
леш еп огпег. 

The machine has a smooth outer shell made of either 
ln cs Jiron or a highly conducting material, such as 
ле о о ено та штит. in order to provide a unitorm boundary 
wid wioOtmam@emcontine the magnetic field Μα Πα Che machine. 
In the case of conductive image shield, circulating image 
ποτ ουσ lo confine the magnetic field: 

Inercia nin highly conductive shield between the 
field winding and armature which rotates with the field 
winding. The purpose of this shell is to intercept 
asynchronous magnetic fields due to load imbalance and 
Space harmonics of the armature field and take large fault 
porques and radial loads under fault HE It 
may also intercept thermal radiation from the room tempera- 
ture armature or a separate thermal radiation shield may be 
used. In any case, the essential purpose of this shell is 
шоп ргогесе tae field windings. Since this shell will not 
be strong enough to take large fault loadings, it should 
possibly be supported by a shell made out of a stronger 
material. In this work, the load carrying structure will 
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be called the shield structure and the shell which provides 
aecual shielding will be called the primary shield. кои 
In kKieures 27and 3, the primary shield is inside the ορ μα. 
structure, but it may very well be outside the shield 
Seructure. Consequences of these configúrations will be 
ersecussed іп а later section. 

The electrical configuration is shown in Figure 4. 
пеге аге а total of six windings, three phase, one field 
and a direct and a quadrature axis damper winding. The 
electrical parameters of the machine are the coefficients 
of the flux-current relationship which can be expressed in 


matrix ρου} 0211. 


^a a Lap Lap Ma “аа аа E 
^p Lap Та tao № Ња (а ip 
Xo E Pa La E Бой eal le 
^" ? ШИ и м TL. Ut. Ὁ P" 
Asa fad) а қа Ha ta ° Та 
Хад Па Пе Ὁ ο. i 
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Figure ! 


ШЕСЕ Еа Ырга ОТ 
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Thi) o may be expressed by: 


= La Lab Lab 
% Lab La Lab 
Ng | Lab bab L. 
Ap — IM cospó M cosp(d-21/3) M cosp(d+21/3) 
А-а М. соврф М. совр(%-21/3) М. соѕр(ф+21/3) 
А-а ШЕ sinpó M. sinp(6-21/3) ШЕ sinp($*27/3) 
M cospó М. соврф Mg sinpó la 
М соѕр(ф-2т/3) ШЕ cosp(d-21/3) ШЫ біпр(%-21/3) ip 
M cosp(d+21/3) М. соѕр(ф+2т/3) М. sinp(ó6*27/3) 1, 
τ Lea $ ір 
Lea Га ° isd 
0 0 La 1sq 


By doing so, it is assumed that all rotor-to-armature 
mutual inductances vary sinusoidally with rotor position. 
This is equivalent to neglecting higher order harmonics. 
It should be emphasized that higher order harmonics in 


phase-to-phase mutual inductances were not neglected in this 


study. 20 
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Field Analysis | 

ror the purposes of this study, a two dimensional 
Пперресепгаршоп is usedi ле ешбесса ог finite length oinn 
епа turn geometry will be included in an approximate fashion. 
In the Stbatent, active section of the machine, currents 
are parallel to the axis of the machine and effective 
armature ana field current densities are related το гета 


Current by 


8 ° L, Nat 
За ‚о тты сы σα (1) 
wae Sao (1-x ) 
µ. Τ᾿ 
RAS E PON (2) 
0 Ре e (1-y ) 


The field and armature winding are assumed to have a 
uniform density between their inner and outer radii, and 
po fill the wedged shape space within a winding angle. 
ENpPressions ког the magnetic fields generated by these 
currents are derived in references 16 and 17, and are 
summarized in Table 1. Wherever double signs (+ or +) 
mppear, the upper sign refers to the laminated iron shield 
case, the lower to the conductive image shield case. Unless 
otherwise specified, all current densities take space factors 
WCO account. 
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TABLE I 


MAGNETIC FIELDS 


ВЕТ Winding Fields 





ГИ 
np?#2 
nd 
27. sin -= sin np(0-49) = np-2 
Ber - r|=— K 
РР We nT(2-np) Res 











2-пр,|2-пр 
- Е E 


ΓΙ; по отрове: 





Ө 4 
Тр sin — зіп 2(0-ф) тыны! Reo 
5 r|-1n ytn(l-y ) RS 


T 





R 
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2-пр, |2-пр 
Б = Fae 


ІЗ 2mp 
2+ 
= η. 





if p=2, п=1 component: 


д 
Jr sin | нге) оя 2(0-9) 1 y. [Reo 4 
= r|-in уір(1-у ) P» 
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Table 1 
Magnetic Fields (continued) 





Ке 
np#2 
——=== | 51п np(8-9) Res np+2 
Пат = — — —r |-2np- (2-np) = 
п ода 
np-2 np-2|R прте T 
+(2+np) —- t(2-np) == = (узу j 
fo S S 








ορ αι component: 


Шы” Rund R EL 
T3 fo fo 
i EN “Πα - o 











np#2 
ES e 
27}, біп |--а--|сов пр(0-%) Ra. Inp+2 

E f 2 Be 

Нер = nm -n?p р -H+(2—np) > 
n odd 
np-2 np-2|R np? + 
+(2+пр) | —| _ +(2-пр) |-Е- πὶ (αρ 
126; 5 5 


ШІ р-2, m=! component: 


сов 2(0-%) | R H R 
г E | 3 +41n ој 
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Table 1 
Magnetic Fields (continued) 





Ray p R. 
ПО fe 
-- 2Ј sin|—5—-|sin one прва, 
pf a m пт(2+пр) p 
LONE mno 
S 

| 2J. sin|—5—-|cos np(8-$) |К, |пр+2 np+2 

Hop” τ nT(2-np IT (1-3 
n саа 


To get Armature phase A fields 


Replace By 
0-0 
Je Ya 
y | X 
wfe Ө ае 
Ho: E 
Προ i 





ШС 


inductance Parameters 

Self and mutual inductances of and between the various 
ΠΡ: επ Calculated Бу integrating flux linkage об the 
eds given in Table 1, over the area of a winding. 
Reference 16 gives such a calculation and the results taken 
from this reference are summarized in Table 2. 

The electrothermal shield has been treated as a 
fourier series of sinusoidally wound damper windings. The 
Justification for this is that each space harmonic field 
will excite currents in this shell that are sinusoidally 
distributed in angle. Further more, each space harmonic is 
Миле сода ο all but one of the ficticious windings: Each 
mailing 15 characterized by a self inductance and mutual 
inductances with field and armature windings. The mutual 
inductance between any two of these damper windings is zero. 
Шис пава in the expressions, Nnp is fictitious number of 
ο αν ог damper winding: 

The expressions given in Table 2 have been left with the 
effective length noted differently for each parameter. 

EE C lengths may not correspond to any physical lengths of 
machine parts, because of the effects of end turns. Further 
more, since end turns of different coils in a machine will 
шеге different lengths, we may expect the coils to have 
differing effective lengths. Certain assumptions utilized 
mien evaluating effective lengths are discussed in section 
"Mot this chapter. 
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. TABLE 2 





INDUCTANCES 
прж2 
16 1 6م‎ N sin! NO ‚Fe 
ШУ аа 5 5 ο. 
n3pm6*, , (n^p*-5) (1-y?) 


R дпр 
k np-2 Eu ШӘ 
vasa [somme s 


но=е. nel component: 





E 3 5 5 5 (пр?) хай е 
п ода п“рпд а р ea -х ) 
2 





Ц пр-2 прта се ο 
(np+2)x ашы а-қ не 





i bD=2, n=l component: 
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М = > M. соз прф 


απ οσα 
where 
ΠΠ ΠΤ... Ν sin ПО ae зла | ute аде пр 
Е TT — m R 
А 2 2 fo 
Du 2 — 2 2 R 
1 О а ΙΤ )(1-x^) а 
апа De 2np 
[1-x97 BP. RP ШЕ В ] 
= np 5 
a 5 Гог прж2 
пр[4- (1р) ] 


gf p=2, n=l component: 


Ето REM ME É 
әз [- тракт я 
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' Ho la n Маре R. 2np 
Lan? —— ] + Б. 1 
бпр 5 
npÉ2 
2u Q N N sin|P*wae cos прф 
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20 GEM | 2 t п-х2798 |520) 
аап 2 | ج‎ E R Мамо - 
пар ا‎ х )(2-пр) ao 
R 2np 
E zs ] 
S 





if p=2, n=1 component: 
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D ] 
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IE 
DN PRO wae ont 
E 5 16 lato at sin | 5 cos |— Cor 
a 3 2 2а 
п ода n pre е Т р -4) (1-х) 
R ПЕ 
Ц пр-2 npt2.2|. ao 
enti orem" 
E B=. Hl component: 
l ἢ N = sum was 
2 alo at 2 (ES “++ s -X 1) 
2 2 8 
T6 P ) 
R 
S 
2 пи ло Ont 
ТО u N sin wae (1-008 (237) 
τ s a o at > = 3 
T6 lex n р 
where 
R 2np 
т 4 221B. аа 
[2-np-lix +(np+2)x Е = ) R ] 
С = 
n 
pU np[4-(np)“] 
for np#2 
Mor p=2, п=1 
dU ee 
C = Ls T οὐ ao| 
52 | aD dO R. 
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Replace 


cos npọ 


COS ROD 


cos прф 





By 


(%ф-2т/3) 


(ф%2т/3) 


(ф-2т/3) 


(ф%2т/3) 


sin npo 


біп пр(%-2т/3) 


біп пр(ф%2т/3) 
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Machine Rating Equations 
ithe volt-ampere rating ora machine IE SIVE Ûy 


Ve 
P = 3V.I, = 3BeI, (gO) (3) 


where Ve» Eolo, аге RMS quantities. 
Тће internal voltage Er is given by 


о) МІ 
Be = — + (4) 
v2 





and 


V 


ο. - 
τα = У1-х 2οος2ψ - x, sinV (see Fig. 5) (5) 





Figure 5 


Voltage-current relationship 
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мћеге X, is synehronous reactance normalized to internal 





voltage 
___ fala, We а "Зар та (6) 
та Er Es 

and 
I Jr fe Reo? (1-y^) 


J 0 в *(1-x?) 
т = а мае ао 


а Nat 


They may be combined to give 























6 0 
_ 21 p+2 р+2_ 2-р_ . | wae| . | wfe 
Р = 5. Wu ln Je J. (1-y Reo R sin|—5 E 5 
ox diii 
zu 
РЁ, (7) 
where Cip IS elven from Tabte 2 
R 2р 
END ερ ao 
(ох у-у 5 Ex R. ] 
Cip Е for p?2 
P[4-p*} 
H [R ] 
EID [ ]j in dp gae ] if pz2 
and X, can be obtained by 
X = d E S sis C 
a В 1 wae|_. wfe uu 
V2 π fo sin 5 sin | y ПЕЛЕ 





E 








where C is given by 


nd 


sinó — 095 u ) 
. : 2 "sn 
n odd n P 
and C snp can БЕ Бато Я {Fo ЕЕЕ 


Transient reactances are calculated with the assumption 
chat: 


a) Field-winding flux is constant 


b) Damper-winding current is zero 


These reactances LEE 
2 
M 
X t = X (1-3 Tr) (8) 
а а 2 (L. Lib Le 
t! = 
Ka Ka 


where Ха is per unit synchronous reactance: 





p ea (9) 


Subtranslent геасёапсеѕ аге саси а сеа а па пе σας 
flux linked by the damper shield is constant. These are 


given ьу20›17 


X " = X." = x (1-2 Ms? | 


2 (=> ee (10) 
а Lab Dal 


Since the shell is symmetric, direct and quadrature axis 


subtransient reactances are the same. 


ша 





ЈЕ Machine Losses | 
There are mainly three different losses in synchronous 


M chines}? 
a) Current-conduction losses: 


ου.” (1-2, 


ΝΕ мае а, ao 
Р = AA (11) 


b) Eddy-current losses: 


2 22 2 > 
30 (1–х ји САВ | а 
p = Wae ao avy ed w (12) 


c ΞΕ 


2 


pcne Ds is the mean-squared magnetic field seen by the 


armature conductors: 
2 1 [ ао 
В gr (13) 
ai 


where Bipot given in Table ТЕ Lhe ото DoD ШГП 
B eG is considered, for a two-pole (one pole pair) iron 
shielded machine: 


θ 
р “wie 
E “τ sinf 2 | % 3 
IM ge o = 
3T 








2 
ап аа ] 
5 


Шел this is substituted into (13): 


6 
ги. ІС Б ЖЫ ΜΠΕ wie (1-у 3)]2 


В ت‎ ааа 
av M 3 
{| (1-х) 2 πο 
=== (1-х oy ee) =+ چک‎ (14) 
SR В 
a al sS 








о) Shield losses: 
If an Image shield is used; "eireulating cumrencts 


will produce losses: 


ћ 
р = _.98 5 5 ШЕР 


where В os is Bg evaluated at shield inside radius, and S 
ls the skin depth of the shield material. 
If an iron shield is used, eddy-currents and hysterisis 


will produce losses: 


2 


E 2 
Popes Eee "a o m (16) 


sh 


where а. is shield dissipation per unit mass of shield 


h 
тега]. The effective lengths will be discussed пп 


Appendix B. 
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ЈЕ 


[Шо Effects 

Although elimination of iron reduces the machine vol- 
ume and weight, it degrades the transient behaviour of 
BENE erPrcondueting machines. Large fault torques and stresses 
will be imposed upon the field winding. In addition to 
fault loading, there will be heat dissipation and field 
current rise, which may drive superconducting field winding 
normal. The protection needed will be provided by a suitable 
shield around the field winding. The expressions for fault | 
loadings under various fault conditions are derived in 
references 1 and 18. For us the maximum values are import- 
ЕШ in order to design a machine whieh can withstand Such 
Bads. 

The maximum radial loads imposed upon the primary 


shield will be under three phase short circuit from load. 


These are given ου 215, 
_ 1 2 2 
ο ο EL (17) 


where H 0 and H 1 are the fields outside and inside the 


shell respectively after the fault: 
Над = -H, [cos (wt ,+¢,-p8 )-cos(o,-p6) J 
- pHalcos (wt ,+6,-p8)-cos($,-p8) J 
-Hof соз (м ,+ф.-рө) (18) 


and 
cos (Wt „+#.-p0) 


of 
R 


where f 








апа Ha is the ammature field at primary shield radius: 











и Qr 72-426 2+p 
b 2. ea" | мае R. ]ρ- -2 [1-х Esas ) ] 
HaT UT ШЕ xem > TA 5 R. Rao | 
(2-p) 
S (19a) 
зу. ер " wae Най 2 
Е E "sin Bou R. [-1n x (1= x be ] (19b) 
d a 


if p=2 


and H is the field generated by field windings at primary 


17 
shield radius: 
9 
: wfe = 
| 29 = sin poer Ra, ΞΡ ур+2 В, гр 
E gp” Re le (1- ур к= = | 20) 
т(2+р) τ LS 


it Should be noted that only the first а а от n D 
have been taken into account. 

When equations (18), (19), and (20) are combined with 
(17) and the maximum stress evaluated when wi „=, 

2 

С + 6. 
o 24H, (1+P)LH, (1+f)+H, glcos” (p & $,.) EL) 

The maximum torque will be applied on the primary 


shield under line-to-line short circuit from load. This 


torque is given py7 >19 
e "2 
Т= T [cos wt (sin wen - sin wt)] (22) 
а | 


and the maximum value will occur at нЕ + т/2 „Гог which 


Tmax 7 153 xcu ο.) 











In above derivations a factor ел”, voltage behind the 


subtransient reactance, was introduced to correct fault 
loads from open circuit to load in an approximate fashion. 


Figure 6 shows a method to calculate е," from which 


E " 
е "трон Их" ова)? + (X," . созује _ (28) 


t 


Т\ХА- η 
Tea 


J ха Та 





Figure 6 
Vector Diagram of Internal, Terminal and Reaction 


Voltages 
EI 








Since, primary shield will not ber strong enough o 
resist fault loadings, it must be supported by a shield 
structure. Stresses induced under radial and torque loads 
are derived in Appendix A. 

Under fault conditions, the subpereconauetine 71ers 
winding may be driven normal, due to the heat dissipation 
uud the field current rise.. Perhaps the most рова MEN 
ο το field current risel® in both de and ac levels. Allow- 
able field current rise will depend on the rated operating 
meine With respect to зџирексопаце овњено тиква тр ercer 
Figure 7 illustrates the allowable field current rise for 
I ven operating point. Ac component of (ome TIGL μμ mh 
rise will be πι. compared to the de component and may 
be neglected. Field current will rise exponentially ana 


ра final value will depend on primary cme ет те Ro o Gb ен 


machine reactances and eircuit breaker opening timer»18, 
A ! 
2er Η τη PUN 
ΑΠ — auc а )] (259) 


where ir and i are the instant and initial values of 


Po 
field current respectively; and, Та" is related to the 
primary shield time constant DE 
xt X ae ae 
а X. Ха _Ха | 51 (26) 
Р | : | О Ко 
лог a thin primary shield, shield time constant : 
l =" 
= = + | — 
en ορ ος 2 R [1 p 1 (27) 


This time constant will also determine the machine swings 
under transient conditions, but this was not included in 
ENTIS study. 
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1G 


Des LEN Constraints 
A superconducting coil is character zed) by Us Sem 


Sample test curve. At a given operating temperature, the 
relation between flux density and current density which can 
be carried without driving superconductor normal, const: tutes 
semen a curve. Any chosen operating point Snould be со Une 
Nor Chis curve. In order to provide stability Taede 
disturbances of electrical, thermal, or mechanical nature 
Which may occur, a material of high electrical and thermal 
Bomauctivity is included in the conductor a) Copocr wie 
usually used for this purpose. High-purity aluminum has not 
Been accepted because of the problem of making a metal- 
fhuecical bond with the superconductor aluncugchn ος Aigner 
eaecectCrical conductivity. Recently a metallurgical pondinas 
pem made** but the brittle Παπ. ο πως 
aluminum as a stabilizer. Figure 8 shows Nb Ti superconductor 
σου sample test data, with the state of the art Pine piven 
meee rerenece 2l. To the τινα, ος CNIS ΠΕ ШШЕ = 
ance may drive superconductor to normal, according to the 
ο τον criteria given in this reference. Although Such 
mene 15 пи1се агр1гагу, 15 аел ες ο ο οι а а 8 
experienced with other magnets. Iti should be pointed out 
tnat the higher flux densities are possible with Nb. on 
conductors, but they are more expensive than Nb Ti 
conductors and cannot be fabricated into complex windings. 

Another major limitation comes from ле saturation 
considerations of the iron shield used to confine magnetic 
fields inside the machine volume. The shield may be locally 
Saturated or the maximum flux density carried by the shield 
may exceed saturation flux density of the iron. These 
considerations will be discussed in detail in Appendix B 
wnen the computer program is described. 

Armature current density will be determined by the 
Basses and cooling method employed. A figure of 3.0 x 199 


amp/m^could be accepted for a nitrogen cooled armature. 
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Superconductor Short Sample Test Data 
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ου πεσε integrity of the machine structure will also 
constitute a problem especially when it is degraded under 
fault conditions. The most critical part of the machine 
WII be the shield struceure Ла” MUSE Domesticos 
arge fault loadings. This will be discussed in Aprenda 
Dd B, in detail. 
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Chapter 2 
RESULTS 


|енсшесекз Used in the Design Study 


For the purpose of this study, calculations were 
performed for a 20,000 HP, two-pole synchronous generator. 
input parameters to the computer program are given in Table 
3. First harmone field in the bore of the field winding 
ı6 Shown in Figure 9 along witn thé ире олсо огы 
sample test curve. Shield structure is assumed to be 
pradsnless steel and low values of working Stresses are 
шлозеп аце to the fact that heat dissipation will degrade 
K erial properties. The 120 degree field winding angle 
mas an advantage of eliminating third harmonic component Ol 
η thus preventing flow of third harmonic Current o Sin 
delta connected machines. In addition to the input para- 
meters given in Table 3, field winding inside radius of 
2.76 inches and outside radius of 5.0 inches were chosen, 
maa the iron shield was assumed vo бе ο ο σπα ος cnd 
ШИИБ (© = 0). 

Although the program could be run several times with 
different field winding radii and input parameters to get 
absolute minimum weight and volume machine, this process 
was omitted here and only a design for a given set of input 


parameters defined in Table 3, was obtained. 
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Barst Harmonic Field in the Bore of the Field Winding 
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ТАВГЕ 3 


DESIGN CONSTANTS 

DESIGN HORSEPOWER- 2050825 

MVA RATING= πο 

POWER FACTOR=1.00 

DESIGN RPM=3600.0 

RATED FIELD CURRENT DENSITY=1.15 10 8 A/M2 

RATED ARMATURE CURRENT DENSITY=3.00 10 6 A/M2 
MAXIMUM FLUX DENSITY IN FIELD WINDING=5.00 TESLAS 
FIELD WINDING ANGLE=120.00 DEG. 

ARMATURE WINDING ANGLE=60.00 DEG. 

ARMATURE SPACE FACTOR=0.27 

ARMATURE FILAMENT DIAMETER=0.040 ШЕ 

YOUNG MODULUS OF SHIELD SUPPORT MATERIAL=29.0 ШЕШЕСІ 
POISSON RATIO OF SHIELD SUPPORT MATERIAL=0. 300 | 
ELECTRICAL CONDUCTIVITY OF SHIELD=6.0 10 7 MHO/M 
DENSITY OF SHIELD SUPPORT MATERIAL=0.290 LBS/IN3 
WORKING STRESS IN TENSION-50000.0 Por 

WORKING STRESS IN COMPRESSION=50000.0 me: 

WORKING STRESS IN SHEAR=30000.0 БО? 

CIRCUIT BREAKER OPENNING TIME=10.0 CYCLES 

RATED FIELD CURRENT=680.0 AMP. 

ALLOWABLE FIELD CURRENT RISE=420.0 AMP. 
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2B 


Effect of Primary Shield Location 

Figure 10 shows the effect of Primary Shield loo on 
where Che induced currents flow Under ondas oa 
Dor different air gap lengths. Caléwlat tons are deme wren 
the same set of input parameters as defined in Table 3. The 
case where the primary shield is inside the shield struct- 
ure, produces less radial stresses than the case where 
primary shield is outside the shield structure, under 
three phase short circuit from load, thus may require dese 
оте ја structure and airgap тееп алал ово үе 
machine volume. Two computer outputs given in Tables 4 
and 5, show this effect. As can be seen, the results are 
шос very different from each other. ОАО л Ее 
where primary shield is inside shield structure gives 
smaller machine volume and length, the direc taco T np D 
оше shield structure to the trapped magnetica ΤΝ,» 
maauee higher heat dissipation thus πεσε... 


refrigeration power required. 
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- TABLE Y 


FIELD WINDING OUTSIDE RADIUS= 5.000 IN. 

FIELD WINDING INSIDE RADIUS= 2.760 In. 

ARMATURE INSIDE RADIUS= 6.451 IN. 

ARMATURE OUTSIDE RADIUS= 9.354 IN. 

SHIELD INSIDE RADIUS= 12.257 ΙΝ. 

SHIELD OUTSIDE RADIUS= 19.887 IN. 

TOTAL LENGTH=58.848 In. 

STRAIGHT SECTION LENGTH- 27.237 IN. 

MACHINE VOLUME= 42.313 FT 3 

WEIGHT OF ACTIVE PARTS= 16418.903 LBS. 

PRIMARY SHIELD IS OUTSIDE THE SHIELD STRUCTURE 

AIR GAP LENGTH= 1.4514 ΤΝ. 

PRIMARY SHIELD OUTSIDE RADIUS= 6.364 IN. 

SHIELD STRUCTURE OUTSIDE RADIUS= 6.364 IN. 

PRIMARY SHIELD THICKNESS= 0.0 In. 

SHIELD STRUCTURE THICKNESS= 1.3410 IN. 

SHIELD DEFLECTION TOWARDS FIELD WINDING UNDER THREE PHASE 
FAULT= 0.0220 In. 

SHIELD DEFLECTION TOWARDS ARMATURE WINDING UNDER THREE 
PHASE FAULT=-0.0205 IN. 

SYNCHRONOUSREACTANCE P.U.=0.2533 

TRANSIENT REACTANCE P.U.=0.1875 

SUBTRANSIENT REACTANCE P.U.=0.1283 

MACHINE LOSSES PERCENT OF MACHINE RATING=0.751 

MAX. RADIAL STRESS UNDER THREE PHASE FAULT= 3482.8 PSI 

MAXIMUM TORQUE UNDER LINE TO LINE SHORT CIRCUIT FROM LOAD 
P.U.= 10.303 
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TABLE 5 


FIELD WINDING OUTSIDE RADIUS= 5.000 IN. 

FIELD WINDING INSIDE RADIUS= 2.760 ΙΝ. 

ARMATURE INSIDE RADIUS- 6.306 IN. 

ARMATURE OUTSIDE RADIUS= 9.282 ІМ. 

SHIELD INSIDE RADIUS= 12.257 ІМ. 

SHIELD OUTSIDE RADIUS- 19.887 IN. 

TOTAL LENGTH= 57.950 IN. 

STRAIGHT SECTION LENGTH= 26.773 IN. 

MACHINE VOLUME= 41.667 FT 3 

WEIGHT OF ACTIVE PARTS=  16169.683 LBS. 

PRIMARY SHIELD IS INSIDE THE SHIELD STRUCTURE 

AIR GAP LENGTH- 1.3063 IN. 

PRIMARY SHIELD OUTSIDE RADIUS= 5.026 IN. 

SHIELD STRUCTURE OUTSIDE RADIUS: 6.2941 IN. 

PRIMARY SHIELD THICKNESS= 0.0 ΙΝ. 

SHIELD STRUCTURE THICKNESS= 1.2680 IN. 

SHIELD DEFLECTION TOWARDS FIELD WINDING UNDER THREE PHASE 
FAULT= 0.0245 ІМ. | 

SHIELD DEFLECTION TOWARDS ARMATURE WINDING UNDER THREE 
PHASE FAULT=-0.0229 IN. 

SYNCHRONOUS REACTANCE P.U.=0.2519 

TRANSIENT REACTANCE P.U.=0.1849 

SUBTRANSIENT REACTANCE P.U.=0.1657 

MACHINE LOSSES PERCENT OF MACHINE RATING=0.749 

MAX. RADIAL STRESS UNDER THREE PHASE FAULT= 3377.9 PSI. 

MAXIMUM TORQUE UNDER LINE TO LINE SHORT CIRCUIT FROM LOAD 
Е. 0.=8. 061 
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2C 


Comparison Between 20,000 HP Conventional Generator and a 


puperconductine Generavor 

The results obtained in previous section are compared 
with a 20,000 HP conventional propulsion generator’, 
Figure 11 shows such a comparison in machine dimensions. 
An additional 15,000 lbs are added to the weight of the 
superconducting generator in order to take account of the 
Table 6 


lists some of the information available about the conven- 


machine casing, bearings excitation system, etc. 
tional generator. Unavailability of additional information 
limits comparison in machine dimensions and weight. As 

баг as refrigeration weight and volume arezeonzernedszen ve 
Will not be anything to say unless rotor 18 Structurally 
and mechanically designed. In addition refrigeration 
system may be integrated to provide for other superconduct- 


mae electric machinery. 


TABLE 6 


20,000 HP Conventional Generator Characteristics 


Maximum Continuous 2 Minute 
Rating Capacity 

KVA 1552530 
RPM 3,600 900 
Volts 4,700 823 
Power Factor 1.0 
Frequency, Hz 60 
Amperes 300% Max. 
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Chapters 
CONCLUSIONS AND RECOMMENDATIONS 


One major conclusion that can” be dean from the Ρος πο” 
meoults is that Süperconducting тасш тел muc Pro 
withstand fault loadings and still have low weight and 
ша оте. 

since most Of the weight comes irom πω μπα r 
meplacement of 1t by a conduetave image али АКО ам ΡΞ. 
to additional weight savings, although image shield may 
mwroduce additional losses. 

Replacement of copper by aluminum as a stabilizer will 
also reduce field winding weight if a suitable fabrication 
technique can be found. 

Conversion of the computer program for many pole 
machines will increase its usefulness and provide informa- 
ton about low speed (100-200RPM) Synebrenom qo On EI 
motors. 

Although a simple technique could be employed which 
@manges input parameters systematically and gives direc Ton 
towards minimum weight or minimum volume machine, it is the 
belief of the author that such technique would require a 
Mat amount of computer time and money  Ὁπιοινποαο 
minimum weight or volume machine may not be best as far as 
тог mechanical design and refrigeration input are 
Esncerned. 

Perhaps, the most questionable assumption is that the 
“ment densities for individual machine parts are uniform, 
especially for the rotor structure. Correct estimation of 
rotor weight can be done only by structural and mechanical 
Мел. Critical speeds of Che госот тау пес Бе acceptable 
Као Chal Case rotor structure OT сиам шом ва be 
шест јеа. Rotor lenght should be ineressed as much as 


possible in order to increase thermal resistance and thereby 


Do 





reduce heat leak into cold region. Саса! греваз инек 
however, limit the rotor length. Inclusion OI above солее 
ations will require longer computation time Оно ia we = 


more useful information about superconduc tire im eNe E 
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Appendix А 
SHIELD DEFLECTION AND STRESSES UNDER FAULT CONDITIONS 


It has been assumed that only the shield structure will 
be effective for carrying fault loads. Figure A-1 shows 
radial loads imposed upon simply supported shield structure 
under three phase fault, with the addition of centrifugal 
loading. It is easier to separate these two loadings and 
superimpose them after each solution is obtained separately. 

Reference 19 gives three differential equations for 


deflections: 





921} 211 vV W 
э 245 962 | ели палат 1 (A-1a) 





Шам u 10 оО оси 
2a 3х96 2 53x? ' a^ 3907 7 97 30 15а2 








2 ڪڪ‎ › = = E = 
É 00 а299 ER | a a: 0 (А-15) 
д р 4 4 
u QU W h Ум 2 W 
es ee = επ. = .ا‎ о + — 
УБ = тэ а 3x* а ар EE 





= 0 (A-1c) 


simply supported boundary conditions imply that, 


ас х = Огапан EE 


These boundary conditions can be satisfied if one chooses 
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TX 


и = cos— 2 А соз пө (А-2а) 
u = sinî £ B sin nê (A-2b) 
L n 
и TX E 
w = sin т 2 C. 6ο5 пө (А-2с) 


It was shown in Chapter 2-F that the radial stresses 
depend on the angular position along the circumference. 
They are given by equation 21. 


g 


A ο а == 


3 А = 5 (1 + cos 28) x sinit 


DEP 


When equations (A-2) and (A-3) are combined with (A-1), 
and the last two terms neglected in equation (A-1b) because 


they will be small for thin shells compared to the other 


озне 
A [an?+(1-v)A*n*J-B (1+0) лтп]+С, [2улт]=0 (дена) 
АОЛ ој та Јев Е3(1-ъ)т 60 n FC РОД Що (A-4b) 


A, Cor3aJ-B, пл [3+ (2-0) 2+121 | +0 [322+ (72+ 2n 2)] 





— пт? 
^ RR (А-4с) 
where А = = 
E 
.h 
o 
απ. 
12(1-v?) 
бо 
end for n=0 and 2, BP, am the result which comes directly 


from equation (A-3) 
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ine coefficients obtained оу solurten sol co (А-4) 
аге too complicated to write explicitly. The results are 


summarized in Table A-l. The stresses can be obtained by, 


00 E 





X ج‎ (е, %уе4-2(Х,%У4)1 (A-5a) 
2 Е 
Ta e [£g tve , -Z (X, *vX..) ] (A-5b) 
И сару (А-5с) 
х 2(1:9Уу хо) 
where 
"ai 
ια 
_ · (А-ба) 
ει. ο, ы = 
МЕ а 90 a са 
EUM l ou 
W 
X= go (А-ба) 
“а NS _ 
Mes 5G ‘a | a се. (А-бе) 
τυ су 
o A (5 x 7 3x36) (A-6f) 


and z is the distance from the mid-surface (positive 
towards the axis of the shell). Equation (A-6) can be 
combined with equation (A-5) to give stresses in axial and 


circumferential directions. 


f 





TABLE A-1 


Solution of equation A-! 


AUS AS UL VEN FT 


(121 (1+v) 
-B, У 
* [eae ms Anz +(1-v)21] 


ΠΠ -21A ] 


L'n*e, 12۸ 1+ 
aD ° [6 1+ (71-3 та +)1- (22۶ ЈИ | Eja: +n? 





mE O МА 


SS ty) -3 (1-0) 02-242 | : ELS 
[2m*+(1-v)4A?7] . [322+n2(m?+422)2)/[2van] | 


B [2 (1+0) Ат ]-А., . [2n?2+(1-v)4r?] 


2. 
ΝΑΤ 
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расе the Shield will rotate ма ти y lq wina ir a 
eentretugal stresses willl be induced та зло зесла ве 
пе есетлопв. Ine centritugal loaqirng ша №5 επ a ccoo - 


loading due to uniform internal pressure given by: 


ow?ah 


g (A-7) 


Reference 19 gives the solution for a circular cylind- 
Focal shell with a uniform internal pressure тек еп 


Senply supported: 


T 


EM Εξ E Dc Um я то : 21 
" ^ - ShD} cos2« 4 cosha« : BinB(x-5) . sinh8(x-7) 
2cos* . cosh = Ше 8 en 19 
cos2= + coshze * “9BBla-3) * ES . 


where AES 
6 = "πα 


апа с = 


г 
Е 


Equation (A-8) is the deflection in vertical direction. 
Eom symmetry, the component мог the dispriseemenc ΠΕ 
ENCcumferential direction vanishes. The expressions for 


Әле strain components then become 


- du _ 
Ex ах (A-9a) 
= ME 5 
ερ = (A-9b) 
19 
and it can be shown that 


Qi BE W z 
--ῳ | (A-10) 
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Equations (A-5), (A-6),(A-7), (A-8), (A-9), (A-10) 
and Table A-1 can be combined to get total deflections 
and stresses along the shell at any point on circumference. 
The component w of the displacement will be maximum at 
midlength (xsL/2) when 9-20 and 0-27/2 denoted by Wi 
respectively. These two deflections will determine the 


and Wo 
required airgap along with shield thickness. If the 
material properties are chosen, it is possible to calculate 
stresses at these points, top and bottom surfaces. Since 
the shear stresses will be zero, these are actually prin- 


cipal stresses. Maximum shear stresses can be obtained by, 





(A-11) 


Under line-to-line shors circuit from 10395, lorz: m e 
torques are applied on the shield. The magnitude of these 
ШЕЗСПЕв 15 given by 3. The shield structure must Бе ουσ ας 
σαι to carry these stresses, Нор the Sheal SCrEeESS ΠΤ. .᾽ 


=> 


mell, according to the maximum shear theory 


CAFE _ Е e 2 1/2 
z max 2 128] Hnc | (A-12) 


this case “e will be assumed to consist of centrifugal 


suresses only, that is, 


2,2 
g. = W£ (AT 
д g 
and ов T 
хе 7 u δαν 
απ ου 


Equations (A-13), (A-14), and (A-12) may be combined to give 


maximum shear stresses. 
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Application of this fault torque may cause buckling of 
shield. Reference. 20 cives criteria forveritical ο ο 


stresses. 
г. 7 0.85 (1/12) (12v?) 9 Pg (n/a )9/ (a 1,17? 
(A-15a) 
for 10(n/aj Сер/а«з(а/ћ) 2 
т = 0.272. (1-v2) 3" g(na) 7? xc 


eI 


for. L/a»3(a/h) 7^? 
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Appendix B 


COMPUTER PROGRAM LISTING AND DESCRIPTION 


Inezrirstzhaermonie лета in che field лр ы 


two pole (p=1) machine is given by (Table 1). 





rfw | 3T 


Zen 3 
kj eo E 


Ө 
ug Jp sin T sin( 6-9) | 
- а 





5 

The maximum flux density in the r direction can be 
found by taking the derivative of В fu with respect to r 
and setting it equal to zero, which gives a simple result 
for a two pole machine 


po Res (4-17) 


and the maximum flux density in the field winding is 


a wfe R NERO 
B > "RO 2 faye e| a) (A-18) 





max 3T De 


This linear relation between b e and Je defines the 


magnet line which was shown in Figs. 8 and 7. If an oper- 


ating point is chosen (given B and Je) and field winding 


max 
inside and outside radiuses are specified, equation (A-18) 
ша be solved for Н, shield inside radius. Although 

fields generated by the armature phases try to oppose the 


а? there 


fields generated by field winding thus reducing Ва 
will be additional field concentration due to armature end 
turns and higher space harmonics. Thus, calculating Das 


by equation (A-18) is a conservative assumption. 
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After the location of” the iron shield is fixed, the 
first thing СО Бе О κος отима ита вјол такоше 
Pig. A-2 shows the flux 1 тикарес through the iron shield. 
The total flux carried by the shield, for a two pole 


machine 


) 


N 


-(т/2 
ФА =] PPro (Rg) ЈА uy Ge (A-19) 


where the first harmonic component of Bafo comes from 


Maple 1 апа 15 given by 


8 
ο ορ O E 2 
в. = —1 = [e (1-y?) “Ы | (2208 





TIO 


When equation (A-20) is evaluated at r = R. and substituted 


into (A-19), integration gives 


8 
š fe 
ф = ду J sin ү HS 
A- 2 





3T R. 
(A-21) 

The limitations on saturation require that 
Pan © Boat η... 1, (А-22) 


When equations (A-21) and (A-22) are combined, the minimum 


EeMelId outside radius, R ШЕ 








OS 
wre 
7. sinj- R. 3 3 
Е = TB T (VEU Ro (A-23) 
sat S 
where Заре ο τας saturation Πχ density of the iron 
shield. P: - 1.5 οι Was used in the program. 
heocaimsavuration of “Ume Iron shield requires thay 
8 | 
Hu, J, sin sie Reo |3 3 
Boro d.) ^ 3T Rg R. (1-47) € Baa (A-24) 
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Figure A-2 


Flux Pattern in Iron Shield 


Algebraic manipulation of equation (A-23) and (A-24) 
gives 
R _— R <R (A-25) 


This is to say that shield thickness must be less than 
or equal to shield inside radius. Since Ro is fixed from 
equation (A-18), calculated nee may not бл this 
Condition. In this case, Program ве је ο ο τοσα 
design and prints a message saying "NO SOLUTION". Same 
message is also printed when equation (A-18) gives Du less 
than Rao: 

The maximum airgap length can be found by assuming 


that armature inside and outside radiuses coincide, that is 


A = R,-t-k 


max fo (A-26) 


where t is the radial distance between shield inside radius 


and armature end turns. In the program, this n ΤΠ. 


ах 
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divided into n steps (n=50 in this program) and in each 
iteration airgap length A is incremented by А а 1 20. 


Given A, armature inside radius Roi 


Ho MS A (A-27) 


сила пат ове ношње ве раат шы i 


DRE Βα τ t)/2.0 (A-28) 
lhe program then proceeds to sum up the infinite 
series given in Table A-2 of Chapter 2 for La - Lap: Then 

synchronous reactance is normalized to internal s s 

and x, as given by equation(6)is calculated. Table A-2 
πιο Che effective епа те ин а о јен = study”. Equation 
7 1s then utilized to give an expression Поссе Еи en 
section length, which results in a second degree equation 
for the straight section length. This осоо Е 
complicated to write explicitly and will be omitted here. 


otal machine length is givenToy 


A. = GE ο (A-29) 

After all machine dimensions are calculated, the 
program then proceeds to calculate ποτ: ΕΘ. σε. 
Synchronous reactance Xa given by(9)swESIDEMEmPMEPOS OIM rc 
Ха" given by(8)can be calculated by summing the series. It 
Should be pointed out that up to IS Бао scree ca. 
he user is free to choose number of harmonies to be Sunmede 

To calculate subtransient reactance ха , the primary 
shield radius has to be specified. Since for a given air- 
gap length, the range of R. is known and this can be done 
easily by incrementing R. in each iteration mT Tor a Urey 
iously calculated airgap length А, R. 15 
(m=50 used in this study) R. = Res D А (A-30) 
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TABLE A-2 
Effective Lengths 


втпасцге епа Сик леи МЕ = ЕК + нат 
Effective Zeus e CEN иге £ _ 

self Iindberamee Ξ ^ S 
Mutual inductance, field-to- Ὁ = 2 
armature m 

Mutual inductance, field-to- Аға = ў 

damper 

Eddy-current loss X ed = 2 
Conduction loss A = 24+ 24: 
Hield мапа еле turn lena AX, = R τ. 
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Then x." cans bemealculared From M 


Next a is to calculate required primary shield 
thickness, to give time constant such that, field current 
rise will not be higher than the specified value (SUB- 
ROUTINE TIC). Since the maximum field current rise will 


be from 25 


тала 
ік а 74 
= Тро Е os | (А-31) 
max 


which depends on ха-ха. . If the maximum rise is less than 


κι 
the specified rise, "i the required shield thickness or 
tme constant is set to zero. This is to say that the 
Primary shield will ре а very CHin ПЕ ис о сав un 
where induced currents flow. сис τ, IC Va: μα 
phat the time constant ol Gli о ет E IP 

Having зрестаттеа регатаву виле ла Попа ви ЕР РАНИ MEN 
stresses To under three-phase fault from load is calculated 
(SUBROUTINE STRES) by( 21), 

The most important section of the program is SUBROUTINE 
SEIEL wnhich designs the shield structurceaccemii me suo, Vie 
Siven working stresses. The designed shield must satisfy 


the following conditions: 


1) Wy + hy ЕЕ W. = A 


where h, is the primary shield thickness, h is the shield 


structure thickness and W4 and М2 are the deflections in 


two directions. 


> nou ы h yo = a if primary shic imagers 
outside the 
shield struct- 
ure 
t t Пур = а -if primary 
shield is inside 
Shield structure 


R 
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3) De < p 
Se s Ys it 
Ss ~ Ne 


where E S 5 are induced stresses in compression, tension 


. W Ф W е 
Suc" sit SIS 
stresses in compression, tension, and shear respectively. 


Б? 


and shear respectively and W are the working 


1) 2 < Weis stresses induced under line- 
to-line short circuit from 
load 

T Σπ.) 
er x 9 mene 


ουσ 15 principal stresses are Calculaved aiver 
W. апа W. are found by equation (A-4). These are at 820 
and 9=1/2, on the top and bottom surface. All of the stres- 
ses mentioned above are at the midlength where maximum 
Berresses occur. Shield lengtn'wascassumcsm s cu PCS 
Section length plus field winding епа а ан ЕСКЕ Е 
each of them by the requirements specified above, the 
progrm proceeds to check on requirement 4. If, in the 
course of these calculations, any one of the requirements 
is not satisfied, computations are repeated by incrementing 
Баета Structure thickness, пп рес ΙΤ: 1S5 fabled 
mich shield structure. Figure A-=3 "snows μας ЛОП to be 
anned by the program. Since the Small асрар lenech 
Ell give minimum welght or volume macnhinessbne solution 
will return to the main program if all four of the require- 
ments are satisfied. There may be other solutions at 


bigger airgap lengths which we will not investigate. 
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Machine volume is calculated by 


and machine weight by 


E. 2 2 2 2 E 2 
[ен E ). S. p, rn (Roe n )2, ptm (Re, Aoi а Pe 
2 2 2 
AR Εν ο Е 
where 3 
р. - armature density = 4500 Kg/m 
D οἳ iron shield density | = 8000 Ke/m> 
Pe - rotor density = 4500 Kg/m? 
P ¿T end turns density = 4000 kg/m? 
Da shield structure density=depends on program 


SIE 
These average densities are used in the program. 


Hor selected values of n, m and increment of hj, the shield 
structure thickness solution may not converge or may 
converge at a longer airgap length, ΙΙ спева гасна 
small. Smaller step sizes will increase computation time 
end give true convergence, So some sort of a comprimise 
would be done. if the solution docs nee еш еркен еп 


program prints a message saying "NO CONVERGENCE", 
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Dou Scanned by tne του ии 
Main program and subroutines are listed in Tables A-3, 
A-4, A-5, and A-6. Table A-6 gives the sample program 


input which is described below. 


Card 1 Format 3F10.2 


Horsepower Rating col. 1-10 
Power Factor со1.11-20 


Design RPM col.21-30 


Card 2 Formar 2F10.2 


Fu be. ning out Side Rmaqinsa col. 1-10 
Iinehes 
παπα Inside ποσα col.11-20 
ае ае = 


TO 





Card 3 Pormat 3217002 


Maximum flux density in field 
winding, Tesla 


koged Cic id cyrrent density 
10 `, ΠΡ ш 


Fleld winding angle, Of 
degrees 


Card 4 Format 4F10.2 


Rated armature current density 
x 10729 пам. 


Armature winding angle, e 
degrees 
Armature space factor 


Armature filament diameter, 
inches 


Gard 5 <- Format PIOR? 


Shield structure material 
young modulus KIO SD 


Shield structure material 
poisson ratio 


Primary shield conductivity 
x L < И mho/m 


Shield structure density, lbs/in? 
Working stresses in tension реши 


Working stresses in compression, 
psi 


Working stresses in shear, psi 
Card 6 Format 3F10.2 


Circuit breaker opening time, 
cycles 


Барса πεια current; атр 


Allowable field current rise,amp 


{| 


col 


col 


cole 


cols 


cols 


CONE 


COLA 


Colm 


Coll 


c o m 
col. 


col 


c oO 


СО 


σοι. 


col. 


11-20 


21-30 


11-20 


21-30 


31-40 


11-20 


21-30 


31-40 
41-50 
51-60 


61-70 


11-20 


21-30 





Card 7 


Format Il, if | is punched үп ево то ΙΗ 
makes calculation for the ease where primary 
shield is under shield structure; and, if 2 is 
punched, calculations are performededcq ко 66 


where primary shield is on the shield structure. 


Format F10.2, distance between shield inside 


radius and armature end turns, inches cod. 11-20 
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